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Abstract 

We give a construction of a family of (weighted) graphs that are pairwise 
cospectral with respect to the normalized Laplacian matrix, or equivalently 
probability transition matrix. This construction can be used to form pairs of 
cospectral graphs with differing number of edges, including situations where one 
graph is a subgraph of the other. The method used to demonstrate cospectrality 
is by showing the characteristic polynomials are equal. 


1 Introduction 


Spectral graph theory studies the relationship between the structure of a graph and 
the eigenvalues of a particular matrix associated with that graph. There are several 
matrices that are commonly studied, each with merits and limitations. These limi¬ 
tations exist because graphs can be constructed which have the same spectrum with 
respect to the matrix and are fundamentally different in some structural aspect. Such 
graphs are called cospectral. 

There are many possible matrices to consider, and the matrix we consider in this 
paper is the normalized Laplacian (see Bi). The rows and columns of this matrix 
are indexed by the vertices, and for a simple graph the matrix is defined as follows: 





if i — j, and vertex i is not isolated; 

if i~j; 

otherwise; 


where di is the degree of vertex i. 
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In this paper we want to look at the more general setting of edge-weighted graphs, 
i.e., there is a symmetric, non-negative weight function, w(i,j) on the edges. The 
degree of a vertex now corresponds to the sum of the weights of the incident edges, 
i.e., di = The normalized Laplacian for weighted graphs is defined in 

the following way: 




( 
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if i = j, and vertex i is not isolated; 
if i~j; 


otherwise. 


(A simple graph corresponds to the case where w(i,j ) € {0,1} for all i,j.) We note 
that when the graph has no isolated vertices, £ can be written as £ = D~ l /‘ 2 (D — 
A)D~ 1//2 , where = w(i,j ) and D is the diagonal degree matrix. Finally, we point 

out that this matrix is connected with the probability transition matrix D~ 1 A of a 
random walk. In particular, two graphs with no isolated vertices are cospectral for £ 
if and only if they are cospectral for D~ l A. 

There has been some interest in the construction of cospectral graphs for the 
normalized Laplacian. Cavers [5] showed that a restricted variation of Godsil-McKay 
switching (see [8]) preserves the spectrum, while Butler and Grout [1] showed that 
gluing in two different special bipartite graphs into some arbitrary graph resulted in 
a pair of cospectral graphs. In both cases, the operation preserved the number of 
edges in the graph. 

On the other hand, it is possible for graphs with differing number of edges to be 
cospectral with respect to the normalized Laplacian. The classic example of this is 
complete bipartite graphs K pq which have spectrum {0, l^" 1 " 9-2 ^, 2} (here the expo¬ 
nent is indicating multiplicity). For example, the (sparse) star K^n-i is cospectral 
with the (dense) regular graph K nn . Until recently, this was the only known con¬ 
struction of cospectral graphs with differing number of edges. Butler and Grout [JJ 
gave some examples of small graphs found by exhaustive computation that differ in 
the number of edges, including some where one graph was a subgraph of the other. 
Butler [2] expanded on this example to form an infinite family and showed how to 
construct many pairs of bipartite graphs which were cospectral. 

In this paper we introduce a new construction of cospectral graphs for the nor¬ 
malized Laplacian which can differ in the number of edges. The basic idea is to form 
a ring of linked modules, and then a similar graph where we interchange the role 
of two of the modules (what we term “toggling”). The resulting pair of graphs are 
cospectral with respect to the normalized Laplacian. An example of this construction 
is shown in Figure [EJ Note that the left graph is a subgraph of the right graph. 

In Section [21 we give a formal description of this family, of toggling, and state the 
main result. In Sectional we show how to compute the characteristic polynomial of 
the normalized Laplacian by using decompositions. We then break the decompositions 
of a graph in our family into those which contain a “long” cycle (see Section [4]) and 
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Figure 1: A pair of cospectral graphs for C related by toggling 


those which do not (see Section [5]), and in particular conclude the characteristic 
polynomials are equal so the graphs must be cospectral. In Section [6] we show how 
to go from weighted graphs to simple graphs which are cospectral with respect to the 
normalized Laplacian. 


2 Construction 

Our family of graphs are formed as a ring composed of three different types of 
(weighted) modules: the path on four vertices, the cycle on four vertices, and the 
edge on two vertices, which we label as P, C, and E, respectively. The modules are 
shown in Figure [2] where we have marked the edge weights using a parameter k where 
k > 0 is for now arbitrary. Each module has special vertices marked with “+” and 
” which can be thought of as poles of a magnet to indicate how consecutive mod¬ 
ules will connect. In particular, the “+” vertex on one module will connect with the 
” vertex on the next module. We will refer to these two special vertices as the 
signed vertices. 


k 2 



k k 

/ 

G , C 

: < 

> , © 


°k+ 1° 


Figure 2: The P, C, and E modules, respectively 


A graph in our family is formed by connecting r modules together in a cycle. In 
particular, such graph can be associated with a word using the letters P, C and E. 
As an example, starting with the top module and reading clockwise, the two graphs 
shown in Figure |T] (where k = 1) have the words PPCCPPPC and CCPPCCCP. Note that 
given a graph in our family there are many possible words, i.e., we can choose any 
module to start and any possible direction. On the other hand, given a word, there 
is a unique graph. 
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Definition 1. Given a word W = l\l-i ■ ■ - l T where G 6 {P, C,E} and r > 3. Then 
G(W) is the graph obtained by connecting the corresponding r modules in cyclic 
order as indicated by the word where consecutive modules connect on the signed 
vertices, and where the final module will connect to the first module. 

We note that the two words we constructed for the graphs in Figure [Hare related 
by interchanging the roles of P and C. This will generalize as follows. 

Definition 2. Given a cyclic word W composed of the letters P, C, and E. Then the 
toggling of W is W T , the word formed by taking W and replacing every P by C and 
every C by P. The occurrences of E are unchanged. 

The motivation for the use of the word “toggling” is to notice that the difference 
between G(W) and G{W T ) is adding or removing the edge on a module which goes 
between the non-signed vertices. In essence, we are switching the states of these 
edges. 

We can now state our main result. 

Theorem 1. For a word W of length at least three using the letters P, C, and E, 
G{W) and G{W T ) are cospectral with respect to the normalized Laplacian. 

We note that if W does not contain the same number of occurrences of P and 
C, then the number of edges in G(W) and G(W T ) will differ and are clearly non¬ 
isomorphic. Among other things, we can construct cospectral simple graphs which 
differ by exactly m edges by setting k — 1 and using a word in P and C where there 
are m more instances of C than of P. There are also some special words W so that 
G(W) is a subgraph of G{W T ). One example of this behavior is W = CC... C and 
W T = PP ... P, though others exist (see Figure CO) • 

3 Computing the characteristic polynomial 

Our approach will involve showing the characteristic polynomials of G(W) and G(W T ) 
are equal. We start by determining how to compute the characteristic polynomial by 
the use of generalized cycle decompositions (see HD- For an n x n matrix M = [rri h] \, 

det(M) = ^ sgn(a)mi )<T (i)m 2l<7 ( 2 ) ■ sgn^WM^)- 

(T^iSn , , (?GSn 

Let Gm denote the digraph which corresponds to M, meaning it has i—>j if and 
only if rriij ^ 0. We can consider a permutation o which contributes a nonzero term 
to det(M). The factors of wm(c) correspond to n edges such that each vertex has 
in-degree and out-degree equal to one, as each vertex will appear as the first and 
second index somewhere in wm(ct). Such a collection of edges is a generalized cycle 
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decomposition of Gm- There are three possible structures in a generalized cycle de¬ 
composition: loops (a directed edge that goes into and out of the same vertex), edges 
(pairs of directed edges i—>j and j—ti), and longer directed cycles. More generally, 
if we think of loops and edges as cycles of length one and two, respectively, then a 
generalized cycle decomposition is a collection of disjoint cycles so that every vertex 
is in exactly one cycle. 

In the case when the matrix M is symmetric, many of these generalized cycle 
decompositions will contribute the same factor to the determinant. For example, 
changing the orientation on a long cycle gives a different decomposition but does not 
change sgn(cr)u;M(c r )- With this in mind we consider decompositions. 

Definition 3. Let G be an undirected (weighted) graph. Then a decomposition , D, 
is a subgraph consisting of disjoint edges and cycles. 


When M is symmetric, we can treat Gm as an undirected graph. Every generalized 
cycle decomposition now corresponds to a unique decomposition, D , by removing 
loops and dropping the orientation on the long cycles. Conversely, if we let s = 
s(D) denote the number of cycles of length at least three in the decomposition D , 
then each decomposition corresponds to a collection of 2 s different generalized cycle 
decompositions. Namely, any vertex not in an edge or a cycle has a loop added, edges 
become cycles of length two, and each of the s cycles of length at least 3 have one of 
two possible orientations chosen. 

If we let e(D) count the number of cycles in the decomposition which have an 
even number of vertices (including edges), and F(D) be the set of isolated edges in 
the decomposition D, then we have the following result. 


Proposition 1. Let G be a weighted graph on n vertices without loops or isolated 
vertices. Then the characteristic polynomial of the normalized Laplacian matrix is 


= y^l)e(J>) 2 «(P)(t - n {i,j}£E(D) w (hj) ri{i,j}eF(£>) w (hj) 


D 


ieV(D) di 


where the sum runs over all decompositions D of the graph G. 

Proof. The characteristic polynomial with respect to the normalized Laplacian can 
be written as 


pit) = det (tl — C) 

= det (tl - D~ 1/2 (D - A)D~ 1/2 ) 

= det ( (t - 1)1 + D~ 1/2 AD~ 1/2 ). 

'--- . - " 

=M 

The graph Gm (ignoring loops) has the same edges and non-edges as G, and so we 
can use decompositions to compute the determinant. 
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In particular, every decomposition of G will relate to 2 S ( D ) generalized cycle de¬ 
compositions. For each such generalized cycle decomposition corresponding to a per¬ 
mutation cr, we have sgn(a) = (—l) e( ' D h We will have n — |K(.D)| loops which each 
contribute (t — 1). The non-loop edges z—will contribute w(i,j) / y/dfdj. Now we 
recall that each vertex in a generalized cycle decomposition has one edge coming in 
and one edge going out, and therefore for each vertex i in I '/(D) we will have y/df 
occurring twice in the denominator giving us the dj. Finally, for cycles of length three 
or greater we only use each edge once in the generalized cycle decomposition, but for 
cycles of length two we use the same edge for both directions and so we use the edge 
twice. □ 


4 Decompositions of G(W) with a long cycle 

Proposition |T] shows that we can determine the characteristic polynomial by looking 
at decompositions of the graph. In this section we will consider the collection of 
decompositions of a graph G(W) which contain a long cycle, i.e., a cycle which passes 
through all of the signed vertices in G(W). We denote the set of these decompositions 
as L. 

Lemma 1. Let W be a word of length r with i occurrences of P and m occurrences 
of C. Then for G{W) we have 


^(_l) e P)2 s (£>)(d - iy^-\V(D)\ np,j}g£:(D) W (L j ) npJ}gF(D) W (b j) 

n 


DeL 


i&V(D) di 


(-1 y-\t - l)2( m +b 
2 T ~ 1 (k + l) m+£ 


Proof. Knowing we have a long cycle yields a lot of information about the decom¬ 
position D in G(W). In particular, for a module of type P or E, the decomposition 
will contain only the edge between the signed vertices. These are shown in Figure [3] 
where edge weights have been removed for clarity. 


o o 

I I 


6—6 o © 

Figure 3: Forced decomposition for P and E, respectively 


For a module of type C the situation is a more interesting as there are three different 
options for the decomposition. Namely, that the long cycle passes only through the 
signed vertices; the long cycle passes through the signed vertices and there is an edge 
between the unsigned vertices; the long cycle passes through all of the vertices. These 
three possibilities are shown in Figure |U 
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o-o 


o 


o 


c> 


o 


6—6 


6—6 


6-6 


Figure 4: Three possible decompositions for C in a long cycled decomposition 


Now suppose that among the m modules of type C that precisely h of them are 
the configuration shown on the left in Figure |4j % of them are the configuration shown 
in the center in Figure 01 and j of them are the configuration shown on the right in 
Figure U The choices of which C modules behave in which way is arbitrary. Summing 
over all the possibilities gives the following. 


— iy^-\y( D )\ n{M}€-E(.£>) W (b •?') npj}gF(D) w (hj) 

n 


D&L 


E 

h+i-\-j=m 


ieV(D) di 


2(—l) T 1 (/c + — 1) 


\it 


( 2 (k + 1)) 


-x 




k A 


(-i)fc 4 

(k(k + i)yj V(fc(fc + i)) 


We have 2 s ^ = 2 because there is only one cycle of length greater than three, 
namely the long cycle which contains all the signed vertices. The e(D) will count 
the number of C modules in the middle configuration and possibly the long cycle 
itself. Regardless of the number of C modules in the configuration on the right, 
the contribution from the long cycle to (—l) e ( D ) will be (—1) T_1 . Consider first the 
contributions of isolated vertices and edge weights of the P and E modules. The 
(k + i y-e-m j g |.j ie we ight 0 f th e edges on the long cycle coming from the modules of 
type E. The (t — 1) 2£ accounts for the isolated vertices from the modules of type P. 
Further, the (2(k + 1)) is the product of the degrees of the signed vertices (each such 
vertex has degree 2 (k + 1) as can be seen by noting that in the modules the signed 
vertices have degree k + 1 and then we identify two such vertices). 

ft remains to account for the portions of the decomposotions formed on the C 
modules which are not the signed vertices. The ( h "-j) = m\/(h\i\j\) is the multino¬ 
mial coefficient for how many ways to choose the different module configurations for 
C, and the final three factors are the contributions from each configuration formed 
by accounting for isolated vertices, edge weights, and the degrees of vertices in the 
decomposition (i.e., a pair of isolated vertices or the product of the edge weights over 
product of degrees). Notice that for the middle configuration, the contribution to 
(—l) e( W h as been appropriately grouped. 

Now we can simplify by pulling out the terms which do not depend on the sum 
and cancelling. For the terms in the sum we can use the multinomial theorem to 
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simplify. Continuing the above computation, we now have 


(_i r -l (t - 1} 2 Y fc 4 fc 4 \ m 

2^~ 1 (k + iy+ m V [k{k + l) 2 ) 2 {k{k + l) 2 )/ 

2- 1 (A; + l) £ + m 1 J 2' r_1 (fc + l) m +t • 

The important thing to note is that the expression in Lemma |T| will be the same 
for W and W T because m +1 is invariant under toggling. 


5 Decompositions of G(W) without a long cycle 

Any cycle in a decomposition with edges in consecutive modules would have to go 
through all of the modules to close up. In particular, if there is not a long cycle in our 
decomposition D, then the decomposition is composed of only edges and CVs which 
lie in individual modules. 

We consider what decompositions can happen in a single module and how decom¬ 
positions in consecutive modules interact. The first task is straightforward to carry 
out, and in Tables mm and [3] we show the possible local decompositions for each 
module. To help facilitate the analysis we have grouped the local decompositions by 
which signed vertices (if any) are used. 


Table 1: Decompositions of the P module 


signed vertex 

local decomposition 

contribution 

neither 

9 9 

& --6 

1 

+ 

LI 

-k 

(t — l) 2 (2k + 2) 

— 

9 c 

Q < 

3 

> 

-k 

(t — l) 2 (2k + 2) 

+/- 

9 9 

6 — 6 « 

3 9 

5 —© 

k 2 - (t- l) 2 
(t — l) 4 (2k + 2) 2 


The next part is to understand the transitions between modules, i.e., how local 
decompositions interact. We have already grouped the local decompositions by which 
of the signed vertices are used. We now note that if signed vertices are used in by 
a local decompositon in one module, it influences which of the signed vertices are 
available for use in the next module. This is indicated by the following transition 



















Table 2: Decompositions of the C module 


signed vertex 

local decomposition 

contribution 

neither 

9—9 9—9 

6—6 6—6 

k 2 

(t — l) 2 (k + l ) 2 

+ 

n 

-k 

2{t — l) 2 (k + l ) 2 

— 

9—9 

6---0 

-k 

2(t — l) 2 (k + l ) 2 

+/- 


-1 


4(f — l) 2 (k + l ) 2 


Table 3: Decompositions of the E module 


used 

local decomposition 

factor 

neither 

Q © 

1 

+ 


0 

— 


0 

+/- 

O © 

-1 

1 

to 


matrix with rows and columns indexed by subsets of the signed vertices: 


Q = 


0 + - +/- 

/ 1 1 1 1 \ 

1111 
10 10 
\ 1 0 1 0 / 


Using Q we can now count the number of ways that we can have decompositions 
use the signed vertices in the modules for G(W). This is done using the transfer 
matrix method (see H), and in particular is equal to the number of closed walks in 
the directed graph corresponding to Q which have the same length as the length of 
the word. We need to go one step further and for every module add the contribution 
of the local decomposition. 

This final part is done by adding in diagonal weight matrices where the diagonal 
entries correspond to the contribution of the decomposition for that particular mod¬ 
ule. These contributions are found by (—1) raised to the number of even cycles (i.e., 
edges or C 4 ’s) times the product of the edge weights used in the local decomposition 
(remembering for an edge to use that edge twice), divided by the product of the de¬ 
grees of any vertex used in the decomposition. The only subtle part is handling the 
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vertices which will not be a part of a decomposition in any module. What we do is 
assume at the beginning that every vertex is isolated and contributes a (t — 1) then 
whenever a vertex becomes a part of the decomposition we divide by (t — 1) to correct 
(the choice of this approach is because signed vertices lie in two modules, hence while 
it might not be in the decomposition of one module it could be in the decomposition 
of the other). When there are several possible decompositions in a given case we 
add them together to form the entry for the weight matrix. The contributions were 
previously listed in the tables and become the diagonal entries of the weight matrices. 
We therefore have the following weight matrices. 


X P 


x c 


x E 


(1 0 

n __ 

V (t-l) 2 (2/fc+2) 

o o 

v o o 


o 

o 


-k 

(i-l) 2 (2fc+2) 

0 


0 \ 

0 

0 

(t-l) 4 (2fc+2) 2 / 


(i-i) 2 (fc+i ) 2 u 

n _ dz _ 

U 2(t—l) 2 (fc+l) 2 

0 0 


V 


0 



(1 

0 

0 

0 

\ 

0 

0 

0 

0 


0 

0 

0 

0 


v° 

0 

0 

-1 

/ 

4(t-l) 2 


0 

0 


—k 

2(t—l) 2 (k-\-l) 2 

o 


0 \ 

0 

0 

-1 , 
4(t—l) 2 (fc+1) 2 / 


So for the graph G{£ i£ 2 • • • £ T ), we have the following: 



(_1 )< D )2< D \t - i)”—W(^)l 


D£L 


n {i,j}eE(D) w (^j) II{ij}eF(D) w (hj) 

n iev(D) di 

= {t- 1) |V(GW)I trac e(QX ei QX h ■ ■ ■ QX £t ). 


( 1 ) 


We now focus on rewriting the trace expression in (PQ). To start we note that 
Q = RSR- 1 where 

/ 3 0 0 0 \ 

0 0 10 
0 0 0 0 ' 

\ 0 0 0 0 / 

Combining this with trac e(AB) = ti&ce^BA) we can conclude 


R = 


/ 1 -1 1 
1 -1 0 
\ 1 0 
U 1 -2 


1 \ 

0 

-1 

0 J 


and S = 


(t - trac e(QX h QXe 2 • • • QX £t ) 

= (t - tmce(RSRr ] X £] RSRr 1 X^ 2 • • • RSR~ 1 X £t ) 
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= (t - l)\ v ^ w »\ trace ((SRr'X^R^SR^X^R) ■ ■ ■ (SR~ 1 X £t R)). 


Because S has two rows of 0’s this simplifies the matrices that we have to deal with. 
In particular we have 


SR-'XpR = ( 

< Yp 

v 0 

SR-'XcR = ( 

( Ye 
v 0 

SR-'XeR = ( 

<Y Z 
. 0 


O 

Z c 

o 

Z E 


, and 


where if we let u t — 1 then 
F P = 


u 1 +32 ku *—s /■ -»-' + u;» 4 — Sfcttf+/.•- — n- 
12 (fc+l) 2 a 4 

Sfcf u 1 + .1 6kn i +2/.'~ U 2 + N» 1 +2/, » 2 — //~ + ■-/ J 
24(fc+l) 2 w 1 


-si-II' 

0(A: + 1 )-// 1 

I ''' J » AM ! — N /,- 2 I f I . : —S/r J — ttf 

12 (A- — 1) -1/ ! 


Yn = 


Vr = 


I(i/7u 2 — :i'2ktr — Hi/.- 2 — 11)»~— 8*fc — J — 8/,: 2 » 2 — J.iiA;u 2 -R8A:~' — Si/ 2 -2k- I 

12(fc+l) 2 u 2 6 (/c+1) 2 m 2 

8fc 2 M 2 + 16fcM 2 —8fc 2 +8M 2 +2fc+l — 4A: 2 m 2 — 8ku 2 +4k 2 — 4m 2 — Ak+1 

24(jfc+l) 2 w 2 


12(fc+l) 2 M 2 


16m 2 -1 —8m 2 —1 


12 m 2 


6m 2 


8m 2 + 1 -4 m 2 + 1 

24m 2 12m 2 


Because we can carry out block matrix multiplication, we note that the resulting 
upper left block will be the product of the upper left blocks and that the resulting 
lower right block will be the all zeroes matrix. This allows us to conclude the following: 


(t - 1) |V(GW)I trac e{QX tl QX t2 • • • QX tr ) = (t - l) |y(G(w))l traced Y 2 • • • Y £t ) 


There is no convenient way to find a simple expression for these decompositions 
as we did for the long cycles. However, it suffices to show that the toggled words will 
produce equivalent results, which is what we now show. 

Lemma 2. Let W = £1 i 2 • • • t T and W T = 7172 ... 7 r . Then 

(t - l) |y(G(,y))l trace(F £l F £2 • • • Y tr ) = (t - i)T( G ( wT ))l trace(Y 71 Y 72 • • • Y 7r ). 


Proof. Both sides are polynomials, and so it suffices to verify that the relationship 
holds for t 7 ^ 0 , 1,2 (i.e., if two polynomials agree at all but three points, they must 
agree everywhere). To show that they are equal, we will make use of the following 
special matrix, 


f 20 u 2 - 2 —32u 2 - 4 \ 

\ 8 m 2 + 1 - 20 m 2 + 2 ) ‘ 


This matrix has the following special properties, which can be verified by carrying 
out matrix multiplication: 
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• UYp = Y C U. 

• UY C = Y P U. 

• UYp = YU. 

These properties are key, in that they indicate we can pass U through one of the Y ,* 
matrices but we need to change the matrix in the same way that we do in the toggling 
operation. 

For t ^ 0,1, 2 we have that U is invertible and so by repeated application of the 
above properties we have 

(t - l) |y(G(H0)l trace(y, i y, 2 • • • Y ir ) = (t - l)\ v (G(w))\ tra c e(UY tl Y t2 ■ • ■ 

= {t- l) |y(G(iy))l trace(F 71 UY^ 2 ■ ■ ■ Y ir U~ x ) 

= (t- l) |y(G(w))l trace(y 7i r 72 C/ • • ■ Y<\ T U~ X ) 

= (t - l) |y(G(iy))l trace(y 71 F 72 • • • UY k \ T U~ X ) 

= {t- l)\ v ^ w ^ trace(y 71 F 72 • • • Y 1t UU~ 1 ) 

= {t- 1)1^^))! trace(y 7i y 72 • • • Y lr ) 

= (t - l)l y ( G ( wT ))l trace(F 7i y 72 • • • Y Jt ), 

where in the last we use that toggling does not change the number of vertices in the 
graph. □ 

Proof of Theorem [0 To show that the graphs G(W) and G{W T ) are cospectral we 
can show that they have the same characteristic polynomial. We use Proposition [Tj 
and consider all the possible decompositions. Lemma [1] shows that the sum over all 
the decompositions which contain a long cycle are equal while Lemma [2] shows that 
the sum over all the decompositions which do not contain a long cycle are also equal. 
Thus the sum over all decompositions is equal, and the theorem is established. □ 

6 Weighted Graphs to Simple Graphs 

We have considered graphs with edge weights in terms of a parameter k as shown 
in Figure [2j By letting k — 1 and restricting to P and C modules we will produce 
cospectral simple graphs. 

Simple graphs can also be obtained by appropriately “blowing up” our graph. 
This works by replacing vertices by independent sets. An edge between u and v 
which has been replaced by r and s vertices respectively then becomes a complete 
bipartite graphs between the two independent sets with all edge weights w(u,v)/rs. 
(Note that r and s are generally chosen so that this new edge weight is 1, i.e., so the 
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new graph is a simple graph.) Similarly several consecutive E edges with weight k + 1 
can become k+1 parallel paths. A discussion on how eigenvalues for the normalized 
Laplacian work for blowups can be found in [2J. In particular, it is known that the 
eigenvalues of the blowups are determined from the eigenvalues of the original graphs 
(which we have shown to be cospectral) and the remaining eigenvalues will come from 
the blowup procedure, which will be the same for both graphs. 

As an example in Figure[5]we start with the cospectral graphs EEEPCC and EEECPP. 
This figure also contains the blowups which result by replacing the unsigned vertices 
in C and P modules with k independent vertices (marked by putting k inside the 
vertex and making the lines bold to represent complete bipartite graphs), the three 
consecutive E edges become k + 1 parallel paths of length three. In particular, the 
resulting blowups are simple graphs which are also cospectral. 

k+1 k+1 





Figure 5: The graphs EEEPCC and EEECPP (above) and their respective blowups (be¬ 
low). 

There are other possibilities. For instance, from the definition of the normalized 
Laplacian we note that the matrix does not change if we scale all edge weights by a 
fixed amount. So we can first scale the edge weights and then perform a blowup. A 
partial example of this is shown in Figure |6] where we consider the graph corresponding 
to ECC. By setting k — 1 and then scaling all edge weights by 2 we get a weighted 
graph which has as a blowup the graph shown on the right in Figure |6j By a similar 
process we could also do the same for EPP to construct a cospectral pair of simple 
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graphs. 



Figure 6: The graph ECC and corresponding blowup, where k — 1 and all edge weights 
have been scaled by 2. 


7 Conclusion 

Many, if not most, approaches to establish cospectrality rely on showing that a small 
perturbation in the graph corresponds to a small, controllable perturbation in the 
eigenvectors and hence eigenvalues are preserved. This was not the case in this 
construction, which is why we considered the characteristic polynomials. Also, while 
we show that the characteristic polynomials are equal, we never explicitly computed 
one. Instead, we showed that the method to determine these polynomials will produce 
the same answer for a pair of cospectral graphs. It would be interesting to find 
additional families where this can occur. 

We have been able to establish a large family of cospectral graphs for the nor¬ 
malized Laplacian (and hence also probability transition matrix) which have unusual 
properties, including cospectral graphs with differing number of edges and graphs 
cospectral with subgraphs. There is a vast amount about the spectrum of the nor¬ 
malized Laplacian that is not well understood. We hope to see more of this area 
explored in future work. 
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